Hole boring in a DT pellet and fast ion ignition with ultra-intense laser pulses 
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Recently achieved high intensities of short laser pulses open new prospects in their application to 
hole boring in inhomogeneous overdense plasmas and for ignition in precompressed DT fusion tar- 
gets. A simple analytical model and numerical simulations demonstrate that pulses with intensities 
exceeding 10 22 W/cm 2 may penetrate deeply into the plasma as a result of efficient ponderomotive 
acceleration of ions in the forward direction. The penetration depth as big as hundreds of microns 
depends on the laser fluence, which has to exceed a few tens of GJ/cm 2 . The fast ions, accelerated 
at the bottom of the channel with an efficiency of more than 20%, show a high directionality and 
may heat the precompressed target core to fusion conditions. 

PACS numbers: PACS: 52.50.Jm, 52.75.Di, 25.20.-x 



A progress in short-pulse lasers nowadays allows one to 
achieve focal intensities larger than 10 20 W/cm 2 , where 
the radiation pressure becomes the dominant effect in 
driving the particle motion [l|, [2|, [3[ . The plasma electrons 
are pushed steadily by this force, and ions are accelerated 
in the strong electrostatic field forming a shocklike struc- 
ture 0, The use of circularly polarized laser light 
might further improve the efficiency of the ponderomo- 
tive ion acceleration avoiding strong electron overheating. 
It allows one to obtain a quasi-monoenergetic ion bunch 
in a homogeneous medium by adjusting the laser pulse 
and plasma parameters 0, 0, M @, 0, & @] • 

In this Letter, we investigate the effect of radiative 
ion acceleration for the first time in inhomogeneous over- 
critical plasmas with densities rising up to more than 
100 n c . The critical electron density n c for laser wave- 
length A = 0.8 /im equals 1.72 x 10 21 cm -3 , which corre- 
sponds to the mass density for a deuterium-tritium (DT) 
plasma of 7.2mg/cm 3 . We show the possibility of hole 
boring in such plasmas driven by laser pulses with inten- 
sities exceeding 10 22 W/cm 2 , which is much deeper and 
more stable than considered before [loj]. Following the 
idea of fast ignition in thermonuclear fusion [ll[ , we pro- 
pose to use fast ions from the outer target layers directly 
irradiated by the laser, to ignite the precompressed fuel 
core. The plasma density profile of the DT target at 
the stagnation moment of the precompression phase was 



obtained in numerical simulations [l2(]. It corresponds 
to the baseline all-DT target of the HiPER project [l3[ 
that is expected to release ~ 10 MJ with the energy gain 
~ 60. The electron density region (1 — 100) n c , where 
the hole boring process is supposed to take place, can be 
roughly approximated by an exponential profile with a 
spatial scale of ~ 20 /im. 

Simulation results described below demonstrate that 



the interaction of an ultra-intense laser pulse with a high- 
density plasma takes the form of a quasi-steady structure 
schematically shown in Fig.[lJ The laser ponderomotive 
potential sweeps all electrons forward, so that the charge 
separation field forms a double layer propagating with 
the velocity Vf, in which ions are accelerated forward. 
The parameters of this double layer (electrostatic shock 
or piston) standing in the range, [#^,0], slowly evolve in 
time due to variations of the plasma density and the laser 
intensity. The laser penetration is bounded by a peak of 
the electron density. Much heavier ions penetrate to the 
x' < region, where they are reflected at x' a and accel- 
erated by the electrostatic field in the forward direction. 
In front of the piston, the plasma is practically neutral 
- ions and electrons reflected from the piston are propa- 
gating in an unperturbed plasma with the same velocity. 
The front shock at x' s = Vft' separates the reflected flow 
and the unperturbed plasma. 

Let us consider the momentum conservation in the 




FIG. 1: Schematic structure of the piston and the electrostatic 
shock maintained by the radiation pressure. The frame is 
moving with the piston velocity. Curves: laser intensity, blue; 
electron density, green; electrostatic field, black; ion density, 
red. 
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frame of reference co-moving with the piston. A Doppler- 
shifted laser beam is coming from the left and reflected 
back from the electron density spike. The momentum 
flux deposited by the photons is 21' / 'c, where the light 
intensity in the piston frame V is related to the in- 
cident laser intensity I by the Lorentz transformation 
V = 1(1- /3f)/(l + fif). This light pressure has to be 
balanced by the pressure of the particles coming from 
the right side with the piston velocity Vf = (3fC (they 
are at rest in the laboratory frame) and will be elasti- 
cally reflected back. The particle momentum flux reads, 
2n / i Vfrrii^fVf = 2riirriiC 2/ y 2 f3 2 ^ where 7/ = (1 — (3 2 )~ 1 ^ 2 
and rii is the upstream ion density in the laboratory 
frame. Here, we neglect the momentum of electrons be- 
cause of their small mass. Then, by equating the photon 
and ion momentum fluxes, one finds the equation 



1 1-/?/ 



niiriiC 3 1 + (3f 



(i) 



for the piston velocity [14]. It can be solved explicitly: 

f = B/(l + B), (2) 

where B = (I /riimiC 3 ) 1 / 2 . In particular, for B <C 1, 
the motion of the piston is non-relativistic, and one finds 
the known formula [H, (3f ~ B. The piston velocity 
characterizes also the reflection coefficient of the laser 
light R = (1 — j3f)/(l + /3f). The piston slows down while 
propagating into denser plasma regions also resulting in a 
decrease in laser absorption: A = 1 — R. The dependence 
of the piston velocity on the plasma density ([2]) is shown 
in Fig.El(a). 

The velocity of ions reflected from the piston in the 
laboratory frame is related to Vf by the Lorentz trans- 
formation Vi = 2/3fc/(l - 
of accelerated ion reads: 



/3j). Accordingly, the energy 



2m iC 2 j 2 f Pf. 



(3) 



In contrast to previous works, let us now consider an in- 
homogeneous plasma with a slowly varying density pro- 
file. Integrating the equation for the piston coordinate 
dxp/dt = Vf over the densities ranging from ni mm to 
nj max , one finds the time needed to push the plasma 
up to the density ni max , X5 = c _1 J drii Lj /3fUi, where 
L(rii) = rii / (drii / dx) is the density scale length. In case 
of an exponential profile with a constant scale length L, 

T b = 2L (m iC /I) 1/2 (n^ 2 ax - n^ 2 in ) + L p /c, (4) 

where L p is the plasma layer thickness. The first term 
in Eq. (|4j) determines the laser pulse duration needed to 
penetrate the chosen plasma layer T p = T^ — Lp/c. Its de- 
pendence on the ion density in the interval (1 — 100) n c is 
plotted in Fig.[2](b). The time of 100 laser periods corre- 
sponds to 270 fs for the laser wavelength 0.8 /im, which is 



a relatively short time compared to the typical lifetime of 
the plasma channel. The latter is about 30 ps long. It can 
be estimated dividing the channel diameter ~ 10 /im by 
the characteristic sound velocity in the ambient plasma, 
< 300 /im/ns, for a temperature of 1 keV. 

Now we can find the energy spectrum and the fluence of 
accelerated ions. Within the time interval dt, the piston 
passes through the distance Vfdt and accelerates dNi = 
noiVfdt ions. These ions are distributed over the energy 
interval de = (dsi/dx) Vfdt. Thus, we find the following 
formula for the ion energy spectrum: 



dNj 
de 



IL 



2mfc 5 /^ 7 6(i +/ ^)' 



(5) 



where (3f is related to the ion energy according to 
Eq. (|3|) . Integrating the ion energy over the density range 
[ni min, Tii max] with a constant scale length L, we find an 
expression for the ion fluence from this plasma region: 
Fi = J Si rii dx. For sufficiently large densities, such that 
(3f <C 1, we find: 



t-, 2/ L Tii max 

Fi ~ In . 

C Tii min 



(6) 



To verify our analytical model and to confirm the the- 
oretical possibility of hole boring in an inhomogeneous 
plasma using superintense laser pulses we performed one- 
dimensional (ID) and 2D simulations with a relativis- 
tic electromagnetic particle-in-cell (PIC) code (l5| . The 
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FIG. 2: (a) Piston velocity as a function of plasma den- 
sity for a circularly polarized laser pulse. Each curve cor- 
responds to the value of dimensionless vector potential a — 
(I /meTicC 3 ) 1 ^ 2 , shown in the legend, (b) Density dependence 
of time T p required for the laser pulse to traverse a plasma 
layer with an exponential density profile with L — 20 /im, 
rii min = n c , and rii max = 100n c . The color of curves corre- 
sponds to the same field amplitudes as in (a), (c) Ion energy 
as a function of plasma density for several values of a. (d) 
Stopping range of ions with energies shown in (c) for the same 
density values. The dashed line represents the maximum areal 
density of the fusion pellet to be heated by ions. 
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code does not include binary collisions between particles, 
but it accounts for the electron radiation losses, which are 
important at high laser intensities [l6j]- We consider the 
interaction of a circularly polarized laser pulse with the 
intensity I = 4 x 10 22 W/cm 2 and wavelength A = 0.8 /xm 
(a = 100) with a deuterium plasma layer. The laser pulse 
has a step-like envelope with a 2A rising edge. At the in- 
stant t = 0, the pulse passes the plasma boundary at 
x = 0. The plasma layer has an exponential profile rais- 
ing from 5 up to 100 n c over a distance of L p = 60 A with 
the density scale length L = 20 A. This exponential layer 
is extended by a uniform plasma layer of a density of 
100 n c to avoid the reflection of accelerated particles. 
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FIG. 3: Interaction of a circularly polarized laser pulse of in- 
tensity 4 x 10 22 W/cm 2 with plasma having an exponential 
profile. Momentum distribution p x (x) for (a) electrons and 
(b) ions at instant t = 30 A/c. (c) Ion energy and (d) momen- 
tum distributions at the final instant T& — 250 A/c. The thin 
continuous line represents the analytical distribution function, 
Eq. ((5}. Simulation parameters are given in the text. 

Figure [3](a) shows the longitudinal momentum distri- 
bution p x (x) for electrons at instant t = 30 A/c. The 
laser driven piston zone at x = 14 A, where the laser 
pulse is reflected, can be easily identified. Most electrons 
are driven by the ponderomotive force of the laser field in 
the forward direction. However, due to the generation of 
a strong charge separation field, part of the cold electrons 
is accelerated backwards. The presence of electrons with 
positive and negative momenta in the piston zone demon- 
strates this fact. The backward accelerated electrons are 
interacting with the incident laser pulse, losing their en- 
ergy by generating the radiation, and eventually reverse 
their motion 16]. For the given simulation parameters, 
~ 1% of laser energy has been transformed to electrons 
and ~10% into high energy photons. The radiation losses 
prevent the electrons from escaping the piston zone and 
thus stabilize the laser propagation. The pulse duration 



T p = 188 A/c needed to penetrate to the density 100 n c 
estimated from Eq. (j4]) agrees well with the simulation. 

Ions driven by the charge separation field are acceler- 
ated and move forward as shown in Fig. [31(b). Because of 
relatively low electron energy, an almost complete neu- 
trality is maintained in front of the piston. In a denser 
plasma, the piston speed decreases, which causes a de- 
crease in ion energy observed in the momentum and en- 
ergy distributions at the final time in Fig. 131(c)- (d). The 
ion energy distribution agrees well with the analytical 
estimate following from Eq. ([5|) [thin line in Fig. [3](c)]. 
The energy fluence carried with ions in the shown energy 
range (50 — 300) MeV amounts to 5.4GJ/cm 2 , which is 
27% in respect to the total laser fluence Fi = 20 GJ/cm 2 . 

Results of 2D PIC simulations are shown in Fig. [H 
They were performed for the same laser and plasma pa- 
rameters as in the ID case and with a flattop transverse 
laser intensity profile with a width of 20A and exponen- 
tial wings. The ion density distribution at two instants 
demonstrates an efficient hole boring in the plasma, a 
clean and a stable channel. The fllamentation process 
is strongly suppressed due to radiation losses, similarly 
to the ID case. The radiation process plays a positive 
role in the laser pulse channeling and ion acceleration, 
as it allows to maintain the electron thermal energy on 
a relatively low level and prevents the electron backward 
motion through the pulse p~6( ] . 
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FIG. 4: Channel formation in plasma and ion acceleration in 
the laser piston regime. Ion density distribution at instants 
(a) 90 A/c and (b) 190 A/c. (c) Angular energy distribution 
and (d) the energy distribution of accelerated ions in the do- 
main \y\ < 10 A. Simulation parameters are the same as in 
Fig. El 

The final angular energy distribution and the energy 
distribution of accelerated ions in the central region 
(I y | < 10 A) are plotted in Figs. H(c) and[4|(d), respec- 
tively. Compared to the ID simulation, the hole boring 
velocity is slightly higher due to the pulse focusing in the 
channel. The major part of ions demonstrates a narrow 
angular distribution with the opening angle less than 6°. 

Let us determine the ion and laser parameters nec- 



essary to ignite the precompressed fuel in the pellet 
core. First, we deduce the ion characteristics. The 
scaling laws [17j define the ignition energy E ig and the 
hot spot radius as functions of the fuel density: 



E ig ~ 18 (po/Pmax) 1 - 85 kJ and r ig ~ 20 (p //w) - 97 /im, 
where po = 300g/cm 3 is the reference density. Then one 
can estimate the energy flux needed to ignite the fuel 
Fig = Eig/Trrfg, which weakly depends on the fuel den- 
sity: 



^^1.4(p max /p r U9 GJ/cm 



(7) 



This energy must be carried by ions with the stopping 
range of the order of the areal mass density in the com- 
pressed fuel pm&xR — 1 — 1.5 g/cm 2 . Calculations of 
the ion stopping range with the Monte-Carlo code TRIM 
[Tit show that it can be interpolated by a power law: 
pli — lO -3 ^* 8 , where U is the ion stopping length. Here 
the ion energy is taken in MeV and the areal density 
pli in g/cm 2 . For the areal density of the compressed 
fuel ~ 1.5 g/cm 2 , the maximum energy of fast ions fully 
stopped in the core is limited by the value ~ 60 MeV. 

The laser and plasma parameters that are needed to 
produce the ions with such an energy can be found from 
Fig. [21(c) and (d). The laser pulse with the amplitude 
a = 70 (J ~ 2 x 10 22 W/cm 2 ) would accelerate ions 
to energies less than 60 MeV in the plasma region with 
densities larger than ni m [ n = 70 n c . Laser pulses with 
smaller amplitudes a < 50 are probably not appropriate 
for stable hole boring, because the long time X5 implies 
a strong sensitivity on 2D effects. The laser pulse may 
lose its propagation stability, if the ponderomotive force 
is not strong enough to evacuate most of the electrons 
and to clean the channel in the underdense plasma jTj|. 
Then, the Raman scattering would destroy the pulse and 
impede its access to the denser plasma layers. 

In estimating the ion fluence, one has to account for 
the divergence of the ion beam. Taking the density of 
compressed core p ma x = 400 g/cm 3 , we estimate the ra- 
dius of ignition spot as 15 pm. Assuming the ion beam 
angular divergence of « 6° taken from Fig.[4](c) and a 
distance of propagation « 50 pm [l2| , one finds that the 
beam radius increases by 5 pm. Then one has to deliver 
the laser energy in the focal spot of 10 pm and to generate 
the ion beam with the fluence Fi ~ 3 GJ/cm 2 . 

Applying Eq. (|6|), we obtain the upper density value 
n>i max — 200 n c of the plasma acceleration layer with an 
exponential density profile L = 20 pm. The ion energy 
spectrum generated in the density range from 70 to 200 n c 
extends from 15 to 60 MeV. The laser pulse duration re- 
quired for the ion acceleration is 0.75 ps, according to 
Eq. dH). Consequently, for the intensity 2 x 10 22 W/cm 2 , 
the laser fluence is 15 GJ/cm 2 . Taking the focal spot 
area 3 x 10 -6 cm 2 we estimate the laser energy of 45 k J. 
This value is of the same order of magnitude as in other 



schemes of fast ignition, although the required power is 
60 PW. Such a high power is supposed to be achieved in 
the second stage of the ELI project [2o| . 

In addition, we have to account for the energy needed 
to assure laser propagation to the acceleration zone 
^imin = 70 n c through the undercritical plasma. Here 
one can use a longer laser pulse with a lower intensity 
rsj 10 21 W/cm 2 . Applying the scaling law from [19[, we 
estimate the pulse duration of 3ps necessary to reach 
the critical density. According to Eq. (|4)) , one needs 4 ps 
more to reach the bottom of acceleration zone. There- 
fore, the hole boring requires the laser fluence of about 
7 GJ/cm 2 , half of the acceleration amount. This value 
can be reduced by appropriate shaping of the laser pulse. 

In conclusion, we have demonstrated analytically and 
with PIC simulations that there is a range of laser pulse 
intensities ~ 10 22 W/cm 2 that could be efficient for hole 
boring, in situ ion acceleration, and fast ignition of pre- 
compressed thermonuclear targets. The advantage of this 
approach is that it does not need any additional arrange- 
ments and can be applied for spherical targets, which are 
well adapted to high repetition rate operation. 
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BLAN07-3- 186728 and by the Region Aquitaine, project 
34293. 
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